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SUMMARY 

Application of a general method for calculation of inviscid 
hypersonic flow fields to the title problem is discussed. It is 6on-' 
eluded that the complications caused by the requirement for sonic flow 
at the rear corner and particularly of the uncertain position of the 
stagnation streamline lead to sufficient difficulties with convergence of 
iterations that a practical procedure is not likely to be found. 


INTRODUCTION 

During the past few years, the numerical analysis of high speed 
inviscid flows has reached a point where numerous codes, reflecting 
a variety of approaches, are available. These vary widely in complex- 
ity and in the precision of their predictions. Among the more simple is 
a method originated by Maslen^ which bears a strong family relation- 
ship to the familiar methods of integral relationships. This method was 
originally intended for analysis of hypersonic axisymmetric flow over 
smooth bodies and it does indeed provide simple and accurate results, 
quite useful for repeated engineering calculations. 

2-7 

Subsequently, a number of extensions have been carried out 
including application to radiating gases and to nonequilibrium chemistry. 
Also, Schneider has presented® an analytic discussion of the basis of 
the method. Another tack^» 10 has been concerned with extension to 
three dimensional flows. In this case the question has not really been 
whether the basic idea is applicable in three dimensions. Rather it has 
been whether the resulting computer code would remain sufficiently 
efficient to be of practical value. 

The present study has been concerned with application of the 
method to a deceptively difficult problem. This is the high speed flow 
over a blunted cone having a large opening angle and at angle of attack. 
The shape is typical of planetary entry vehicles, notably the Viking 
Aeroshell. The geometry is simple; the problem, of course, is that 
the velocity near the surface remains subsonic for interesting cone 
angles, reaching sonic values only at the rear corner. 
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ANALYSIS 


A. General Considerations 

Inasmuch as the main analytic structure of the method has 
been given elsewhere9, 10 , the development will not be repeated 
here in detail, although Appendix A shows the final system of equations 
used. A flow diagram and identification of the program elements is 
given in Appendix B. 

The basis of the analysis is this. Consider the equations of 
motion in a shock oriented system, and observe that the shock layer is 
thin. The variation of pressure across the layer (controlled by the 
lateral momentum equation) Is smooth and generally moderate as com- 
pared to, say, the density variation. The latter is related to the 
entropy variation behind the curved shock and generally falls rapidly 
toward the body, so that the mass flow in the shock layer is concentrated 
near the shock. On the basis of these arguments, it is arbitrarily 
assumed that the variation in pressure across the layer is quadratic in 
the stream function > i. e. , in the mass flow - and can be written 
explicitly in terms of the pressure at the shock and its normal gradient. 
Thus 


P = P 


shock 


+ Pj (*- 4-shock) + P 2 



- \} shock 
2 ^ shock 
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The sum (P^ + P^ ) is the normal gradient; the split is based on arguments 

for simple cases. Then for a given shock, the entropy is known on a 
streamline and the remaining thermodynamic properties are directly 
calculable as is the speed for an isoenergetic condition. Thus for axi- 
symmetric or plane flow, there remains only to invert the definition of 
the stream function to obtain the physical coordinates. For unsymmetric 
flow, an ordinary differential equation (eq. All, Appendix A) must be 
solved first to determine the streamline direction (laterally). All this is 
quite straightforward and good results have been obtained in a number of 
cases. 


The complications arise almost entirely from the real requirement 
that the body, not the shock, be specified so that one is led to an iterative 
process. For axi symmetry, this iteration is not difficult and can be per- 
formed rapidly. For some three-dimensional situations the same statement 
applies even though the requirement to solve a differential equation hurts. 
The main question is the number of iterations required. 

For the present case of a blunt cone of large opening angle, there 
are two severe complications each of which separately very much in- 
creases the number of iterations required. One is the occurrence of sonic 
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flow at the rear corner and the other deals with the position of the 
stagnation point. We examine these in turn. ' 

B. Sonic Corner , . 

Consider first the logic used for, say, a hemisphere cylinder. 

The procedure^ has been to postulate the shock curvature as a quadratic 
in distance along it and to iterate this form to satisfy simultaneously the 
nose position and curvature and the body position and slope at a point (A) 
downstream of sonic flow. The actual iteration proceeds easily. Subse- 
quently the solution continues in a marching fashion with a very straight- 
forward iterative process. So far, no problem. The result is quite 
insensitive to the position of point (A). On the other hand, for the present 
geometry, the sonic point is at the rear corner so that the above logic 
would force one to do the whole flow with a quadratic shock curvature. 

This is much too restrictive. The alternate procedure is to choose A at 
some forward position and try to go ahead. What happens is familiar 
from the method of integral relations. The solution typically diverges 
before reaching to the rear end. By altering the position of the match 
point (A) in an iterative way, the solution yielding the desired sonic 
corner value can be found. Such a process works well but obviously 
increases the computing time. When the ratio of specific heats (y) drops, 
convergence becomes slower and the sonic behavior becomes a much more 
localized phenomenon. Parenthetically, it may be noted that McDonald^ 
has given an alternate procedure which should work well, though some of 
his results for low y are not correct. 

In three dimensions, it should be apparent that this sonic con- 
dition gets messy, and may indeed be a mixed condition! 1 , with the flow 
supersonic on a portion of the rear edge. One may well require that the 
(no\^ line of the match point (corresponding to point A discussed earlier) 
might have to vary in downstream location depending on the angular 
position. Such an event boggles my mind. 

C. Stagnation Region 

For axisymmetric flowthe position and shape of the stagnation 
streamline are trivially given. For an unsymmetric case with a rounded 
geometry giving supersonic flow quickly, the position and shape are not 
hard to find in the spirit of this method. Particularly for high Mach 
number, the shock is loosely parallel to the body and it has been assumed 
that is exactly so and thus that the stagnation streamline is straight and 
normal to body and shock. Even though this may not be strictly true 
the assumption usually appears to introduce no important error. 

Now consider the case of interest here. Figure 1 shows a 
spherically blunted cone of 60° half angle, at 20° angle of attack. At 
point B, the body is normal to the flow direction. If the corresponding 
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shock is normal at the associated point (K), we are led to a shape 
like, say CKD. 

Unfortunately, it is quite reasonable that the shock be more like 
EFG for which the normal point (F) is far removed from B. Further- 
more there is no longer much reason to expect a simple stagnation line 
(FH?). The complications grow apace . 


D. Effects on Calculation 

So long as the match line (A) can be at a uniform axial distance 
downstream, the sonic corner condition causes no problem beyond that 
described for the axisymmetric case- Figure 2 shows calculated re- 
sults from the present program run for a zero angle of attack case and 
compares these with those from an axisymmetric program and a result 
given by South^ . The iteration proceeded without trouble. 

When the same configuration is at angle of attack, as in Figure 1, 
disaster strikes. The program has no objection to computing the flow 
behind either of the shocks sketched. What it does get confused at is 
the idea of an iterative process converging on what is now, from the 
shock viewpoint, an entirely arbitrary body shape. Efforts to date to 
conduct such an iteration have met with no success. 

Numerous calculations were attempted within the framework of 
shocks like CKD(fig. 1); that is, shocks whose normal point (K) lies 
opposite the corresponding body normal point (B ) so that the stagnation 
line (KB) is parallel to the freestream. The results showed clearly 
that convergence of the iterative process cannot occur within this limited 
family of shocks. Next, calculations were extended by simply moving 
the normal point (F, fig. 1) an arbitrary distance relative to point B. 

The results showed considerable improvement over the previous cases. 
However convergence was not obtained, nor did a reasonable iteration 
logic appear. Recall the logic used^ in the absence of this complication. 
We assume a shock whose curvature (CS) is given by 

CS = A(0 ) + zB(0) + z 2 D(0) 


where the most general form of A is 

A(0) = - 5 

1 + b Cos^e 


Z and 0 are cylindrical coordinates, with origin at the shock normal 
point ( K , fig. 1). Now a, b are found by matching to the body 
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position, slope and curvature at the point opposite K. Then B(0), 

D(0) are found by matching position and slope at a Z = constant surface 
downstream. Even this logic requires some compromise to permit 
reasonably rapid convergence. Now we add a displacement of the shock 
to F so that the body opposite no longer need be normal. 

While it seems that, in principle, the desired convergence can 
be obtained, the process appears so tedious as to be foolish -- remember 
that the whole object of this study was to provide a reasonably rapid 
means of computing the flow field. 


CONCLUSIONS 

The assumptions made to simplify the system of equations of 
motion do indeed lead to a flow model capable of computing the flow 
behind a quite general smooth shock in three dimensions. However, 
the practical problem of finding the shock for a specified body leads to 
an iterative process which can greatly increase the computing time. 

For the geometry studied in the present work, the very large cone 
angle leads, at appreciable angle of attack, to a non-negligible displace- 
ment of the stagnation point. This unfortunate circumstance introduces 
another dimension to the iterative process making it so extensive as to 
preclude its practical use. 
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APPENDIX A 


The Euler equations of motion are to be solved for a 
compressible, isoenergetic flow. Consider a set of coordinates (z,rj, 
8) related to the flow geometry (Fig. Al). We base them on the 
shock where z, 0 are the axial and azimuthal coordinates of a point on 
the shock while r\ is the (normal) distance from shock to field point. 
These coordinates are related to cylindrical ones (x, r, w ) by 

x = z +r) Sin T f 

+r)^Cos^T - Zr\ r Cost Cos X (1) 

6 S 

r Sin (co- 0 ) = r| Cos t SinX 

where r = r (z, 0) is the equation of the shock and v , X , and t are 

3 3 

angles (Fig. -Al) defined by 


Tan v 


Tan X 


Tan t 


8rs 

9z 

I JUjL 

r 86 


Tan v Cos X 


Now introduce three mutually perpendicular velocity components, 
U, V Cos £, Y Sin £ where U is normal to the shock and £ is otherwise 
arbitrarily chosen. Then 

U = Uj Sin T " v 2 ^ os T C° s ( w “ 0+X) 4- w^Cos t Sin (a; - 0 + X) 

V Cos | = Uj CosT + v^Sinr Cos(w-0+X) - w^Sinr Sin (u-0 +X) (2) 

V Sin £ = -v ^ Sin (w-0+X) - w^Cos (<a>- 0 +\) 
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We define a pair of stream functions, i|i and <t > , such that the 
velocity vector U is 

pU = V4> x (3) 

and let 4, cr , i{j be the final independent variables related to the physical 
coordinates by 

t = z 

cr — 0 (4) 

i 

^ (n> z > e ) 

Then Eq. (3) yields 

L (*) = 0 


- 6 

T cr 

pV71 > = e(A Cos £ + B Sin £) 
where the operator L is 


(5) 


L = ~ P? s . § ,. + |) JL — + (A Cos i + B Sin £ ) (6) 

r 9 cr 3L, 

where A, B, D, and E are geometric factors given by 
A=[r s +rt Cost Cos X (X 0 - 1 )] / r Cos \ 



T 

e 


D - Tan T Tan \ + Cost 

z 

E = (1 + nr Cos t)/Cos t 

z 7 



The entropy (S) equation is 


L(S) = 0 (7) 

while the isoenergetic equation relates the entropy (H) to Y, neglecting 
the small normal velocity (u, eq. 2) by 


1 

h + 2 V = constant 

The approximate normal momentum equation is, on 
integrating, 


( 8 ) 


P<*, o-) = P s (^,cr) + P 1 (4.o-){^-^ s )-P 2 (^ «r) 



where 


P 


1 


and 




pV 


2 

Cos t CosX[Cos £ L (r)-Sin | Cos t L 



(10) 


| n^ p -£* U^[ Cost -Cos t Cos \(1 g)/r-Sin TSin X Tg/r]| 


The remaining momentum equation describes the turning 
(£, Eqs. (2) ) of the streamlines. One has 

dSinl = l jstn jL - - T - a - p - T gin \ V Cos J Cos \ I 
d £, pv 2 j ^ r *ay r Cos t 


+ Sin £[~ ( ‘-'OsXTan4(Xg-l) _ | Cos £ CosX Tan t 
L r J r Cos t 


( 11 ) 


Equations (4) and (7) - (11) plus a state equation form a 
determinate system. The one differential operation is L which is the 
ordinary derivative along a streamline. After solving, one returns 
to physical space, finding q by Eq. (5). 
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APPENDIX B 


The computer program consists of a main program plus a 
number of subroutines, of which the important ones are: 


CROS (I) Tests the position of the body from TRAN 

against the actual geometric coordinates 
and provides an iteration scheme. 

EXTRAP Extrapolates the solution one step down 

stream when the results in CROS converge. 

INTEGG (I) A third order Runge-Kutta routine to inte- 

grate the equation for the streamline. 

NONAX Main logic program controlling the sequence 

of use of the subroutines, 

SHOKQ (I) Computes quantities at the shock for 

specified radius, slope and curvature 
distribution, 

START Provides initial values at the nose as well 

as starting estimates of the flow near it. 

TEST Routine to test (and provide an iteration to 

improve) the starting solution, 

TRAN (I) Transfer the results of INTEGG from the 

<t , 4* coordinates to z, 0, rj ones and 
integrates to find t) . 

The general logic (shown in the flow diagram) is to make an 
initial estimate of the shock slope and position (START), use that to 
find shock values for 5 stations (SHOKQ), integrate the streamline 
equations (INTEGG) and transform to physical space (TRAN) and test 
the result (TEST, CROS) to determine an acceptable nose curvature 
distribution. Then one goes through INTEGG, and TRAN to station 8, 
uses CROS (8) to test the body position and slope at station 8. Once 
that converges the solution is extrapolated (EXTRAP) and one goes 
thru SHOKQ, INTEGG, TRAN and CROS, iterating as necessary. 



FLOW DIAGRAM 
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(downstream portion) 
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